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Abstract. Wc evaluate the result of the recent pioneer- 
ing numerical simulations in Spitkovsky 2006 on the spin- 
down of an oblique relativistic magnetic dipolc rotator. 
Our discussion is based on our experience from two ideal- 
ized cases, that of an aligned dipole rotator, and that of 
an oblique split-monopole rotator. We conclude that the 
issue of electromagnetic pulsar spindown may not have 
been resolved yet. 



1. Introduction 

We believe that pulsars are spinning down magnetized 
neutron stars with non-aligned rotation and magnetic axes 
(oblique rotators). We also believe that pulsars lose rota- 
tional energy through electromagnetic torques due to the 
establishment of large scale electric currents in their mag- 
netospheres. Unfortunately, an exact expression for the 
electromagnetic (Poynting) flux energy loss L as a func- 
tion of the magnetic inclination angle 9 remains still elu- 
sive. Most people are content with the approximation that 
pulsars spin down at the same rate as 90° vacuum dipole 
rotators, namely that 



L{9) «L vac (90°) = 



B 2 r 6 n 4 
6c 3 



(1) 



where, B, r, fl are the polar value of the magnetic field, the 
radius of the neutron star, and the star's angular velocity 
respectively. Knowing the exact electromagnetic energy 
loss rate, observations of the pulse period P and period 
derivative P would allow us to accurately infer the value 
of B. Moreover, the dependence of L on the magnetic in- 
clination angle 9 has very important implications for the 
evolution and distribution of pulsars in the P — P diagram 
(Contopoulos & Spitkovsky 2006). 

In a real pulsar, the vacuum approximation leading to 
eq.[T]is obviously not valid, because without the establish- 
ment of magnetospheric electric currents a magnetic rota- 
tor cannot generate pulses. The study of the oblique rela- 
tivistic magnetic rotator is, therefore, a formidable MHD 
problem as described clearly in Spitkovsky 2006: "Mod- 
eling of the structure of the highly magnetized magnc- 



tospheres of neutron stars requires solving for the self- 
consistent behavior of plasma in strong fields, where field 
energy can dominate the energy of the plasma. This is 
difficult to do with the standard numerical methods of 
MHD which are forced to evolve plasma inertial terms 
even when they are small compared to the field terms. In 
these cases it is possible to reformulate the problem and 
instead of solving for the plasma dynamics in strong fields, 
solve for the dynamic of fields in the presence of conduct- 
ing plasma. This is the approach of force-free electrody- 
namics" . In his recent pioneering work, Spitkovsky 2006 
managed to evolve numerically an oblique dipolar mag- 
nctosphere for about 1.2 turns of the star, and claimed 
that the solution very quickly settles to a constant elec- 
tromagnetic energy flux which depends on the magnetic 
inclination angle as 

B 2 r 6 ft 4 

£s P itkovsky(60 = (1 + sin 2 9) . (2) 

As we said, this result is very important, and while waiting 
for independent numerical confirmation, we thought that 
we should be able to reproduce its main elements from 
first principles. We will, therefore, consider two idealized 
cases where we know the electromagnetic energy loss rate 
analytically. 

2. The axisymmetric magnetic dipole rotator 

Contopoulos, Kazanas & Fendt 1999 (hereafter CKF) first 
obtained the structure of the axisymmetric pulsar magne- 
tosphere with a dipole stellar magnetic field. They showed 
that, at large distances from the central star, the asymp- 
totic structure is that of a magnetic split monopole, i.e. a 
certain amount of initially dipolar magnetic flux 

irBr 3 fl 

* opcn = 1.23 — (3) 

stretches out radially to infinity in one hemisphere, and 
returns to the star in the other hemisphere. Electromag- 
netic spindown is due to the establishment of a poloidal 
electric current distribution 
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which flows to large distances along open field lines and 
returns to the star through an equatorial current sheet 
that joins at the light cylinder with the separatrix be- 
tween open and closed field lines. CKF obtained a first 
approximation for ^opon and self consistently, by re- 
quiring that the solution be continuous and smooth at 
the light cylinder. The CKF solution has since been con- 
firmed, improved and generalized in Gruzinov 2005a, Con- 
topoulos 2005, Timokhin 2005, Komissarov 2005, McKin- 
nery 2005 and Spitkovsky 2006. 

The electromagnetic spindown luminosity is thus ob- 
tained as 

L C KF = — / J(#)dtf » 3 « ^-^open 

(5) 

(e.g. Beskin 1997). This result may be easily understood 
as follows. As we said, the structure of the axisymmetric 
dipolc rotator magnctosphcrc approaches asymptotically 
that of an axisymmetric split-monopole (Michel 1991), 
and therefore, the spindown luminosity depends only on 
the amount of open field lines $ open . 

3. The oblique magnetic split-monopole rotator 

Bogovalov 1999 showed that eq. [5] is also valid for an 
oblique split-monopolole rotator. As long as current sheet 
discontinuities are present to guarantee magnetic flux con- 
servation, properties of cold MHD plasma flows do not 
depend on the direction of the magnetic field. In fact, all 
properties obtained for the axisymmetric split-monopole 
rotator are the same for the oblique rotator as well. In 
particular, rotational losses of the oblique split-monopole 
rotator are independent of the inclination angle. They de- 
pend only on the square of the amount of open field lines 
'I' open (obviously, the sign of ^open is irrelevant). 

We may now derive an important conclusion about the 
oblique dipole rotator which too, as we argued, becomes 
asymptotically split-monopole- like: rotational losses of the 
oblique rotator depend indirectly on the inclination angle 
only through the amount of open field lines ^open as 

iesthnateW « ^*opc„W ■ (6) 

4. The oblique magnetic dipole rotator 

We have now come full circle. We may estimate the elec- 
tromagnetic spindown luminosity of an oblique rotator 
through an estimate of the amount of open magnetic flux 
as a function of the inclination angle ^ pon{9). A crude 
estimate may be obtained if we calculate the amount of 
magnetic flux that crosses the light cylinder distance in 
the case of an inclined magnetostatic dipole at the origin, 
and rescale our result by the CKF factor 1.23 (note that 
Gruzinov 2005b used a similar argument in his estimate 



of the spindown luminosity). This calculation yields 

*o P cn(0) ~ 1.23^^(1 - 0.2 sin 2 9) , (7) 

which results in the following crude estimate of the spin- 
down luminosiity, 

R2 6q4 

A»timato(0) (1 - 0.4 sin 2 6) . (8) 

In particular, we see that £(90°) < L(0°), contrary to 
what both Gruzinov 2005b and Spitkovsky 2006 obtain. 
Note that Spitkovsky's numerical simulation result agrees 
well with Gruzinov's estimate, and yet the latter mentions 
that he is not even sure about the sign of the 9 term that 
he obtains. 

We do not claim here that our above estimate has the 
same weight as the result of a detailed numerical sim- 
ulation. Interestingly enough, the recent discovery of an 
"on/off" pulsar (Kramer et al. 2006) offers some further 
insight from observations. The pulsar Kramer at al. dis- 
covered is thought to be a 90° oblique rotator with two 
states of emission. During the "off" state, no pulses are 
seen, and the neutron star is indirectly inferred to spin 
down at a certain rate L Q g. During the "on" state, the 
pulsar is observed to spin down at a rate L on — 1.5L Q g. It 
is very tempting to associate the "off" state with a vac- 
uum state where the spindown rate is analytically known 
(eq. [T]) . The observations then imply that 

L(90°) = L on = ^L vac (90°) = £ckf = L(0°) . (9) 

In other words, a 90° oblique MHD rotator is observed to 
spin down at the same rate as an aligned rotator. This 
result is more consistent with our crude estimate in eq. [8] 
than with eq. [2] 

We thus conclude that the discrepancy between the 
two estimates is an indication that the issue of electro- 
magnetic pulsar spindown may not have been resolved yet, 
and that further investigation is needed, both analytical 
and numerical. 
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